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SYMBOLS 
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" skin friction coefficient 

■ Mach number 

■ pressure 

- distance normal to centerline 

■ radius of the duct 

■ total stress tensor 

• time averaged velocity in primary flow direction 

■ time averaged velocity normal to centerline 

* distance parallel to centerline 

- R-r 

■ boundary layer thickness 

■ molecular viscosity 

- y/t 

■ time averagea density 

■ shear stress 


Subscripts 

e ■ boundary layer edge condition 

00 ■ free stream condition 

Superscript 

<( )’> = time averaged fluctuation value 
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Za the study of supersonic turbulent boundary layer flow the turbulent 

shear stress distribution hss always been of great Importance and interest. 

The direct nessurenent of the turbulent shear stress is, however, quits difficult. 

A natural alternative is to conpute the shear fron experimental seen flow data 

by numerically integrating the aonsntun equation. 8uch computations have been 

1 2 

performed in recent studies by Bushnell and Morris , Horstman sad Owen , and 
8turek 5 . This note describes results obtained by a computational procedure 
which differs from those previously reported in that integrated mass and woman- 
turn flux profiles and differentials of these integral quantities are used la the 
computations so that local evaluation of the streamwiss velocity gradient is not 
necessary. Tbs computed results are compared with measured shear stress data 
obtained by using hot wire anemometer and laser veloclmeter techniques in racaat 
studies by Rose end Johnson*' 5 . The measurements of Roes and Johnson were made 
upstream and downstream of an adiabatic unstpa rated interaction of an oblique 
shock wave with the turbulent boundary layer on the flat wall of a two-dimen- 
slonal, M m - 2.9 wind tunnel. The shock wave was generated by a 7* wedge. 

The turbulence data obtained from the two independent systems of measurement were 

in reasonably good agreement, Indicating that the data should be reliable. The 

. - ■ ■ - 

computational procedure developed here is easy to use and the computed results 

show reasonably good overall agreement with those obtained by direct meaaurement. 

7 . - j 

As would be expected for any method of computing sheir stress from mean flow 
data, the computed values of shear stress are quite sensitive to small differences 
in mean flow profi3.es and to simplifying assumptions wMch may be made In de- 
veloping the relationships to be used in the computations. The effect of some 
of these differences on >ompuced shear stress distributions is discussed. 
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BASIC EQUATIONS AND BOUNDARY CONDITIONS 

The time-averaged equations for the conservation of mass and momentum 
for steady compressible turbulent boundary layer flow In an axisymmetric 
channel are, respectively, 

5“ (Pu) + '-P*u’> + 7 gf (rpv) + 7 g7 ( r < P , v*>) = 0 (1) 

and 
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vith representing the viscous stress. 

If we assume that |v <p’u'>| « |p <u’v* >[ and <p’u’>I << | (pu) | 
ana transform to an x-y coordinate system, the continuity and momentum equations 
may be combined and integrated in a direction normal to the surface to yield 
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T » P ■j™ •” P < tt’V 1 >, 

Iquatlon (5) becomes applicable to two-dimensional flow as 1 «*• ». 

The normal stress, T^, which appears in equation (5), is not known fro* 

«iaan profile data* However, computations which have been made in this study 
show that its effect is small. In the results shown, the streamwise gradient 
of T has been neglected. 

knowledge Is also required of the static pressure distribution in the 
boundary layer. In many studies of supersonic boundary layer flow no attempt 
is made to measure the static pressure variation normal to the wall, even though, 

for some adverse pressure gradient flows, the variation may, in fact, be rather 

% 

large. In most instances the static pressure at the boundary layer edge may be 
determined with confidence. If this is done, the normal pressure variation may 
then be represented in approximate fashion by assuming a linear distribution 
between the wall static pressure and the pressure at the boundary layer edge. 
Results ere shown here for both a linear static pressure distribution and a 
constant static pressure. 

Examination of equation (5) shows that an accurate value of boundary layer 
growth rate is very important for the calculation of the shear stress. However, 
precise determination of the boundary layer thickness from experimental mean data 
is difficult. It is even wore difficult to evaluate the boundary layer growth 
rata accurately. This problem may be avoided by using the condition that, the 
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sheer stress diminishes to zero at the boundary layer edge and solving 
equation (5) for dS/dx: 
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In equation (6) a linear pressure variation normal to the wall has been assumed. 

If the streamwlse gradient of;; the normal stress is ignored and if the 
flow is assumed to be locally similar, equation (5) reduces to the express* n 
used by Bushnell and Morris’*' in their computation of shear stress. 

Before solving for the shear stress from equation (5), it is also necessary 
to know the coefficient of skin friction. Thin may b« obtained by using the wall- 
wake velocity profile proposed by Sun and Childs”, the method of least-squares 
may be used to fit the wall-wake profile to the experimental mean velocity pro- 
files to provide values of and to provide a smoothed representation of the 

y 

mean velocity distribution. 

RFSULTS 


In carrying out the cewj.r.t^fLrms for the t'.’.ov downstream of the '.Stock wove 
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boundary layer Interactions, the departure from local similarity has been 
taken into account. For purposes of comparison, however, computations based 
on local similarity have also been made. The effect of static pressure variation 
normal to the wall has also been considered for the downstream stations by assum- 
ing a linear variation in pressure. In these computations the static pressure 
at the boundary la> edge has been computed from the free stream total pressure 
and pitot pressure, with appropriate allowance made for the loss in total pressure 
across the shock system. The velocity and density profiles needed for the compu- 
tations were obtained from mean flow pitot profiles with appropriate allowance 
made for static pressure variation across the boundary layer, and under the 
assumption of constant to total temperature across the boundary layer. As was 
mentioned in the previous section, the mean velocity profiles may be smoothed 
by using a least squares fit of the wall-wake velocity profile to the experimental 
profiles. Computations have been made for both smoothed and unsmoothed profiles. 

Figure 1 shows shear stress distributions computed for an upstream (x = 5.375cm) 
station in Johnson’s and Rose's investigation, along with their shear stress data 
from the hot-wire anemometer and the laser velocimeter measurements. The computed 
shear values which are shown have been obtained under th assumption of local 
similarity. As is shown* the calculated results agree quite well over much of 
the boundary layer, with the measured results obtained with the laser velocimeter. 

The differences between the calculated results and the hot-wire results are 
greater. For both the hot-wire and laser velocimeter measurements the peak values 
of shear stress are seen to occur substantially farther from the wall than is 
observed for the calculated distributions. Also shown in the figure are values of 
the wall shear stress as determined by a least squares fit of the wall-wake profile 
to the mean data and as measured by a Preston tube. The agreement between the two 
shear stress values is good. 



8 


For purposes of comparison, a wall-wake profile proposed earlier by Maise 
7 

and McDonald h^s been used to represent the velocity profiles obtained by Rose 
and Johnson. As is shown, the shear values computed from the Malse-McDonald 
profile are substantially lower than those obtained with the unsmoothed data or 
the wall-wake representation described in Reference 6. 

Results of the computations for the data of Rose and Johnson for a station 
downstream of the shock wave-boundary layer interaction (x = 9.375 cm) are shown 
in Figure 2, along with their measured shear stress distributions. Several sets 
of computations have been made from the mean data, including one, for purposes of 
comparison, in which local similarity has been assumed. Results have been obtained 
for both smoothed and unsmoothed profiles and for both linear and constant static 
pressure distribution across the boundary layer. The velocity and density profiles 
used in computing the shear stress are shown in Figure 3, along with the wall -wake 
representations of the velocity profiles. 

As Is shown in Figure 2, the measured shear stress values obtained by the 
two experimental techniques are in quite good agreement. These in turn agree 
reasonably well with the computed values obtained by using smoothed nonsimilar 
profiles and the assumption of a linear pressure variation across the boundary 
layer. There is considerable differencebetween the results for the smoothed 
and unsmoothed profiles, with the peak shear stress value for the unsmoothed 
profile occuring lower in the boundary layer. Near the wall the difference 
apparently was due primarily to inaccuracy in the numerical integration process 
for the unsmoothed data. Much smaller step sizes could be used with the smoothed 
data. The results obtained from the wall-wake profiles are valid, of course, only 
if the profiles provide accurate representations of the actual velocity dis- 
tributions. In the outer part of the boundary layer the wall-wake profile pr>.u *d- 
an excellent representation of the data at one station but did not fit ro '.-31 .-f 



the other. This would account for some of the difference in computed shear stress 
for smoothed and unsmoothed profiles in that region. 

The computed sheer values are quite sensitive to the assumptions regarding 
static pressure. The computed values of static pressure at the boundary layer 
edge at the first and second measuring stations differed froe those at the wall 
by 3.0 and 1.2 percent, respectively. As is shown, the computed values of shear 
stress with no consideration given to the pressure difference across the boundary 
layer were substantially lower than those obtained when the pressure difference 
was considered. 

As is also shown in Figure 2, the results obtained under the assumption of 
local similarity are markedly different from those determined when similarity is 
not assumed. The peak shear stress levels computed assuming local similarity are 
less than half the values computed when similarity is not assumed. Furthermore, 
the shapes of the shear stress distribution curves are quite different. These 
results occur even though the differences between profiles at two closely spaced 
stations are small. The density and velocity profiles in Figure 3 were obtained 
at streamwise stations located approximately one boundary layer thickness apart. 
Although the profiles at successive stations appear at first glance to be quite 
similar, striking differences are found for the computed snear stress distributions. 

COKCLtfSlMiS 

An alternative meshm of computing shear distribution from experi- 

mental mean profile dare in compressible totor., 1 - r>c 'uand.-ry layer f.H*w has been 
developed. The met !«>.•-] in -.Afferent from tni-. r • i .•.»*!>• reported in that 
Integra tea mass anc rc** *. ‘M c ! j :/. profiles am ■ * . ~ ..vne t. sis of these integral 
quantities ar«_ used ^ fupuca' Jons so that .« evaluation of the t ; *eaF< 
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turbulence data for two-dimensional p.iabatic boundary layer flow In the 
regions upstream and downstream of an oblique shock wave interaction. Thu 
computed results are quite sensitive to the accuracy of the numerical inte- 
grations required in the computational procedure and to the mean property 
distributions in the boundary layer. The assumption of Jocal similarity may 
cause large errors in computed shear stress values for flows subjected to 
pressure gradients, even though adjacent profiles of the mean properties may 
appear, on first examination, to be quite similar. The computed shear stress 
levels are quite sensitive to the static pressure distribution normal to the 
wall. Thus, if reliable shear stress distributions are to be obtained from 
mean profile data, the static pressure distribution must be known rather 
accurately. The effect of the streamwise gradient of the normal stress on the 
computed results is small and apparently may be neglected. 
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Figure 1« 


Turbulent shear stress distribution upstream 
of a shock wave-boundary layer Interact 4 on* 
two-dimensional tunnel [4] 


Figure 2. Turbulent shear stress distribution 

downstream of a shock wave-boundary layer 
Interaction, two-dimensional tunnel [5] 


Figure 3 


Density and velocity profiles downstream 

of a shock wave-boundary layer Interaction (5] 
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□ mean data, linear static pressure distribution, 
Eq. (5) 
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